Competition between charge and spin order in the t — U — V extended Hubbard model 

on the triangular lattice. 
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Several new classes of compounds can be modeled in first approximation by electrons on the 
triangular lattice that interact through on-site repulsion U as well as nearest-neighbor repulsion V. 
This extended Hubbard model on a triangular lattice has been studied mostly in the strong coupling 
limit for only a few types of instabilities. Using the extended two-particle self consistent approach 
(ETPSC), that is valid at weak to intermediate coupling, we present an unbiased study of the 
density and interaction dependent crossover diagram for spin and charge density wave instabilities 
of the normal state at arbitrary wave vector. When U dominates over V and electron filling is large, 
instabilities are chiefly in the spin sector and are controlled mostly by Fermi surface properties. 
Increasing V eventually leads to charge instabilities. In the latter case, it is mostly the wave vector 
dependence of the vertex that determines the wave vector of the instability rather than Fermi surface 
properties. At small filling, non-trivial instabilities appear only beyond the weak coupling limit. 
There again, charge density wave instabilities are favored over a wide range of dopings by large V at 
wave vectors corresponding to y^S) x y^3) superlattice in real space. Commensurate fillings do not 
play a special role for this instability. Increasing U leads to competition with ferromagnetism. At 
negative values of U or V, neglecting superconducting fluctuations, one finds that charge instabilities 
are favored. In general, the crossover diagram presents a rich variety of instabilities. We also 
show that thermal charge-density wave fluctuations in the renormalized classical regime can open a 
pseudogap in the single-particle spectral weight, just as spin or superconducting fluctuations. 

PACS numbers: 71.10.Fd, 05.30.Fk, 71.10.-w 



I. INTRODUCTION 

One of the outstanding problems in quantum many- 
body physics is to understand quasi-two dimensional sys- 
tems where both electron-electron interaction and geo- 
metric frustration are important^. The triangular lat- 
tice is the prime example where the geometry frustrates 
near-neighbor ant i- alignment of the spins that naturally 
tends to occur in the presence of short-range electron- 
electron interaction. Studying models of interacting elec- 
trons on such a lattice is thus certainly of fundamental 
interest, but it is also strongly motivated by the discov- 
ery of new materials. Prime examples of these materials 
are organic bis-(ethylenedithio) (BEDT) Cu2(CN)3 lay- 
ered compounds^, triangular lattice antiferromagnets of 
the CuCr02 family"^, and transition- metal oxide materi- 
als like Naa;Co02 and Nai_a;Ti02. The layered cobal- 
tates have drawn much attention because of their uncon- 
ventional properties. Sodium cobaltate shows an unusu- 
ally strong thcrmopower^ at doping a: « 2/3 that can be 
suppressed drastically by applying an in-plan magnetic 
field^. The observation of Curie- Weiss behavior in the 
magnetic susceptibility while resistivity displays metallic 
behavior is another puzzle^. The system also becomes su- 
perconductor when it is diluted by wate r^i^d" . Various 
types of charge- and spin orders also have been found in 
the system for various dopingsiiii^ii^. 

Na2:Co02 consist of two-dimensional C0O2 layers sep- 
arated by insulating Na^^ layers. The C0O2 layers have 
Co atoms at the center of oxygen octahedra forming a 
2D triangular lattice. The band structure calculation 



performed by Singbi^, revealed details of splitting of the 
bands in Co atoms. With help of this calculation 
and also of NMR experiments^, one can find a rough 
estimate of hopping and exchange constants that would 
enter a two-dimensional Hubbard or t — J model for this 
system. However, the modeling is complicated by the 
fact that band structure calculations lead to hole pock- 
ets that are not observed experimentally, a question that 
is still debated by several groups using, for example, the 
Gutzwiller approximation^^, the local density approxi- 
mation plus Hubbard^^ U and dynamical-mean field the- 
Qj, yi7,i8^i9 _ addition, the effect of long-range Coulomb 
interaction from the sodium leads to modifications to the 
simplest Hubbard Hamiltonian for the cobaltatesS2i2i. 

In this paper, we do not address the question of de- 
tailed modeling of the cobaltates or of other triangular 
lattice systems. Instead, we note that since several types 
of spin and charge density waves are observed in these 
materials, it is quite likely that first-neighbor repulsion 
V, and not only on-site repulsion U, must be taken into 
account. U by itself tends to favor spin-density waves. 
We thus just focus on the simplest extended t — U — V 
one-band model Hubbard model on the triangular lattice 
and ask a few general questions: What types of phases 
are typical in different doping ranges, what type of in- 
teraction favors them, and should one expect pseudogap 
effects. 

Previous theoretical and numerical works have ob- 
tained phase diagrams for the triangular lattice in the 
presence of competing interactions. There are, for ex- 
ample, variational Monte-Carlo calculationsSSiM for the 
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extended Hubbard model. That work focused mostly on 
the presence of the Charge density wave (CDW) at filling 
n = 2/3 and RVB superconductivity at n = 1/3. Slave 
boson methods were used for the t — V and t — J mod- 
gjg24|25 |.Q study CDW, ferromagnetism and also RVB 
superconductivity. Series expansion methods and cluster 
mean field theor}*2^ have also investigated CDW, Neel 
order, ferromagnetic order, dimer order and phase sepa- 
ration in a t — J — V model. We will comment further 
on some of these calculations in the context of our own 
results. 

The results of this paper are obtained with the re- 
cently developed Extended two-particle self-consistent 
approach22i2^(ETPSC) that is valid from weak to in- 
termediate coupling. This method has been bench- 
marked against Quantum Monte-Carlo (QMC) simula- 
tions (QMC) for the extended Hubbard model on a 
square lattice. The approach satisfies conservation laws 
and the Mermin- Wagner theorem stating that no con- 
tinuous symmetry can be broken at finite temperature 
in two dimensions. More traditional methods, such as 
the Random phase approximation, do not satisfy this 
requirement. With ETPSC, quantum renormalization 
of interactions (Kanamori-Briiuckner screening) is taken 
into account. Instability towards zero-temperature long- 
range order is signaled at finite temperature by crossover 
to the renormalized-classical regime where the correlation 
length grows exponentially. The wave-vector of the in- 
stability is determined self-consistently by the approach 



and all wave vectors are in principle allowed. No a priori 
selection is necessary. 

Within ETPSC we can also compute the self-energy 
and other related quantities, such as the spectral weight 
that is measured in photoemission experiments^. For 
the Hubbard model, it has been shown with the Two- 
particle self-consistent approach (TPSC) that a pseudo- 
gap can appear as precursor induced by either antifer- 
romagnetio^l or superconducting fluctuations^!^. The 
former^^ has been observed experimentally in electron- 
doped high-temperature superconductors'^'*. Our results 
demonstrate that CDW fluctuations can also induce a 
pseudogap. This is a relevant question experimentally 
given that CDW induced pseudogaps are observed and 
sometimes even show similarities with observations in 
high-temperature superconductors'^^. 

In the following we first introduce the model and the 
ETPSC methodology. We next present our numerical re- 
sults, discussing various physical effects in terms of the 
spin and charge structure factors. We display phase dia- 
grams that help understand how microscopic parameters 
favor various phases. The CDW induced pseudogap and 
its effect on the Fermi surface are discussed before we 
present an overview and a conclusion. 

II. MODEL AND METHOD 

We write the extended Hubbard Hamiltonian in the 
following form. 
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where Cio- (c;^) are annihilation (creation) operators for 
electrons of spin a at site i of a triangular lattice, riia 
is the density operator, and t is the hopping matrix el- 
ement. The quantities U and V are the on-site and 
nearest-neighbor interactions respectively and fi is the 
chemical potential. 

For the Hubbard model {V = 0), TPSC is a very re- 
liable approach up to intermediate coupling limit. The 
functional derivative method is particularly convenient to 
obtain the TPSC equations'^-. This is the method that 
was used to generalize TPSC to the extended Hubbard 
model^l^. 



The equations that need to be solved are the following. 
The charge and spin response functions take the form 



where 
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are the charge and spin vertex functions and x°(q, ujn) is 
the free response frinction (non-interacting susceptibihty) 
given by 



x°(q,wn) 



rfp /"(p+§)-/°(p-§) 

BZ ^ ~ £p+q/2 + £p-q/2 



(5) 



with 



Eq = —2t[cos{qxa) + 2cos{qxa/2) cos(qj^ V3a/2)] (6) 

the non-interacting dispersion relation and 7(q) — 
—eq/2t. In the above formula i> is the volume of the 
BriUouin zone (BZ), /0(q) = 1/[1 + exp((e, - ^io)/T)] is 
the Fermi function and /io is the non-interacting chemi- 
cal potential. The pair correlation functions are related 
to the static structure factors by 



S'cc(i"i) = 1 + - / — [S'cc(q) - 1] exp(iq • ri), 
n Jbz V 

gss{n) = - [ —[SssicO - 1] exp(iq • ri), 
nJBZ V 



(7) 



(8) 



where S'cc,ss(q) — Saa{c\) ± ^^^(q) are the charge and 
spin component of the static structure factor. The spin 
resolved static structure factor is defined by Saa'is^) — 
(nCT(q)no-'(Q)) ^"^^ ^<T(q) is the Fourier transform of 
n\a. The quantities gccio.) and gss(a) entering the vertices 
Eq.Q are simply the pair correlation functions at the 
first-neighbor distance. 

Self-consistency is established by connecting the static 
structure factors to the response functions by the 
fluctuation-dissipation theorem 
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s(q,'^n), 



(9) 



where u;„ = 2m:T are Bosonic Matsubara frequency. 
Substituting the expression for the susceptibilities Eqs.([2] 



and[3]) and the corresponding vertices Eq. ([4]) on the right- 
hand side, one can use the result to obtain the pair cor- 
relation functions gss and gcc entering the vertices using 
their relations Eqs.([7|and[5]) to the structure factors. As- 
suming that the functional derivatives of the pair corre- 
lation functions arc known, as discussed below, we need 
only three equations to determine the pair correlation 
functions entering the vertices. This is because the Pauli 
principle imposes that ga-,cr{Q) — 0. The equation that is 
dropped out is that for gcdO) = 0. This procedure and 
its impact on the Pauli principle is discussed in detail in 
Ref. |2l. 

Functional derivatives of the pair correlation functions 
with respect to density and magnetization enter the spin 
and charge vertices. The functional derivatives are ob- 
tained from the following equations: 



Sgs{l,2) 

5n{l) 
Sn{l) 

Sgsih2) 

Sm{l) 



- [l-5cc(l,2)], 

■ = 2[l-5u(l,l)], 
= [l-5cc(l,2)]. 



(10) 



These equations are strictly valid only when particle-hole 
symmetry is satisfied. On the square lattice, it has been 
checked by comparisons with QMC calculations that the 
results are satisfactory even in the absence of the full 
particle-hole symmetrj«2^. Apparently, particle-hole sym- 
metry due to linearization of the dispersion relations near 
the Fermi surface suffices. We will make this assumption 
for the triangular lattice where strict particle-hole sym- 
metry is not satisfied. This is justified a posteriori by our 
results. Those that can be checked against variational 
QMC, for example, are in excellent agreement. 

Finally, the self-energy needed to address the pseudo- 
gap problem is obtained following Ref. '271 
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where ojn — (2n + 1)ttT is the Fermionic Matsubara fre- 
quency and is the Bosonic one. We can also find the 
spectral function A{q,uj) — ~^G{q,uj)/2TT. The above 
formula does not assume a Migdal theorem since one of 
the vertices is renormalized. However, it takes into ac- 
count only the longitudinal fluctuations. Transverse fluc- 
tuations could be accounted for following a generalization 



of the steps in Ref. [s^- Since the pseudogap appears only 
when fluctuations are large, the longitudinal case suffices 
to establish the qualitative results. 

Finally, the interacting chemical potential is obtained 
from 

n = TY,f — v4(q,c.„). (12) 



FIG. 1: (Color online) The non-interacting Fermi surfaces at 
fillings, starting from the center, n = 0.5, 1, 1.25 and 1.5. 



III. RESPONSE FUNCTIONS 

The "phase diagrams" that we present in the section 
that follows this one, are determined from the growth of 
the the spin and charge response functions as tempera- 
ture decreases. When the interaction is local, U, the wave 
vector of the instability is determined entirely from the 
non-interacting susceptibility, in other words from nest- 
ing properties of the Fermi surface. The introduction 
of near-neighbor repulsion V changes this since it intro- 
duces a wave vector dependence to the vertices. In order 
to disentangle the various effects, we present the non- 
interacting susceptibilities in the first subsection and the 
results with interactions in the second subsection. These 
numerical results are obtained from Eqs. pi fTU)) . 



A. Non-interacting susceptibility 

The non- interacting susceptibility Eq. [5] is determined 
mostly by the shape of the non-interacting Fermi surface 
that is in turn determined by the dispersion relation given 
in Eq. ([5]). In Fig. [T] we present the Fermi surface for in- 
creasing values of the density, n = 0.5, 1, 1.25 and 1.5 
respectively. The first Brillouin zone is plotted as a solid 
line. It is important to notice that the Fermi surface 
corresponding to n — 1.5 (long dash) touches the first 
Brillouin zone boundary and has long parallel segments 
that lead to near nesting. We will show in more details 
that this causes a strong peak in the non-interacting re- 
sponse function. We also draw two wave- vectors that are 
often found for the most important charge or spin density 
waves in the system. At these wave-vectors, the charge 
or spin response functions often have a strong peak. The 
real-space modulations corresponding to these the wave- 
vectors are depicted in Fig. [21 

The non-interacting response function Eq. [5] is drawn 
in Fig. [3] for different values of densities n — 
0.75, 1, 1.25 1.5 and 1.75 at T = 0.2. The largest re- 
sponse is for n = 1.5. While one might have expected 
that parallel segments of the Fermi surface would have 



FIG. 2: (Color online) Real-space structure for two types of 
order: CDWl on the left panel and CDW2 on the right panel 
are related respectively to wave- vectors Qi and Q2 in Fig. [T] 

lead to a peak at a single dominant wave vector, it seems 
that the frustration imposes a less pronounced maximum. 
However, the height of the maximum at that density in- 
creases rapidly with decreasing temperature. The posi- 
tion of all the peaks changes only slightly with temper- 
ature. The height of the peaks for the smaller values 
of density does not change drastically with decreasing 
temperature. That fact in addition to quantum renor- 
malization of the interactions are the main reasons for 
the absence of any instability at low density up to inter- 
mediate coupling. 

There is a deep minimum near the K point at higher 
values of the density. This is the main reason for ab- 
sence of commensurate spin density wave (SDW). The 
non-interacting response function has a peak at the com- 
mensurate wave-vector Q2 (K point) at lower value of 
the density but, as we just mentioned, this peak does 
not grow enough to produce any sort of order including 
SDW up to intermediate coupling. That is not the case 
in the strong coupling limit but that is out of reach of 
our approach in the density regime where Mott physics 
is dominant. 

It is quite remarkable that the free response function 
shows a strong peak at the origin for n = 1.75, a signature 
for ferromagnetism at nearby densities. It seems that, 
as we will see, frustration on the triangular lattice favors 
ferromagnetism at intermediate coupling, contrary to the 
square lattice case. 



B. Interacting response functions 

In the presence of both types of interactions U and V, 
the response functions are strongly modified. Consider 
typical values of the interaction, [/ — 4 and V = 1.5. 
Using the same color code as in the figure for the non- 
interacting case, we show in Fig. |4]the spin, Eq.([3]), and 
in Fig. O the charge, Eq.([2]), response functions at T = 
0.4 for the same densities as in the non-interacting case 
Fig. El 

In the ordinary random phase approximation, the spin 
response function in influenced only by the interaction U 
and the maxima are at the same location as the in non- 
interacting case. In the present approach however, the 
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FIG. 3: (Color online) The free response function for different 
at n = 0.75, 1, 1.25 1.5 and 1.75 for T = 0.2 



nearest-neighbor interaction V also influences the spin 
response, introducing a wave vector dependent vertex. 
Hence, some of the maxima of the spin response function 
in Fig. m are not at the same wave vector as in the non- 
interacting case. Nevertheless, the differences are much 
smaller than for the charge response function appearing 
in Fig. [5l fn the latter case, the position of the maxi- 
mum is near point K (wave vector Q2) for all densities, 
in other words the charge response is dominated by the 
wave vector dependent vertex introduced by V . 

One can summarize the results for the location of the 
maximum spin response function as follows, fn the range 
n < f .5 the tendency is towards an incommensurate spin- 
density wave (fSDW) while the response is maximum 
near zero wave vector (ferromagnetism) just above this 
density. Generally the height of the peaks increases when 
V is reduced, hence nearest-neighbor repulsion does not 
favor spin order. For n = 0.75, the maximum is near 
point K corresponding to the same lattice structure as 
the CDW2 depicted in Fig. [2] except that one should re- 
place the big or small points with up and down spins. 
This resembles the spin structure that would arise with 
ferrimagnetic order. 

For the charge response in Fig. [SI tendency towards 
CDW2 order (K point) is robust for these values of U 
and V . The tendency is strongest at low values of the 
density because of a weaker effect of the frustration that 
leads to a dip in the non-interacting response function in 
Fig.[3l The effect of frustration is very strong for densities 
very close to n = f .5. But nevertheless the height of the 
peak for densities around n = 1.5 grows dramatically 
with decreasing temperature. 

We verify the dominant effects of U and V discussed 
above, this time by fixing the filling at n = 1.5 and chang- 
ing the interactions. Fig.[n]and Fig. [7]show, respectively, 
the effect on the charge and spin response functions at 
T = 0.4. 



r M K r 



FIG. 4: (Color online) The spin response function at fixed 
U = A, V = 1.^, T = GA and difi'erent value of the densities 
n = 0.75, f , 1.25, 1.5 and 1.75. 
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FIG. 5: (Color online) The charge response function at (7 = 4, 
V — 1.5 and T = 0.4 for different value of the densities. 

For the charge response function in Fig. [SJ we imply 
[/ = 4 for those curves where the value of U is not writ- 
ten. A simple comparison of the charge response function 
with the corresponding non-interacting susceptibility for 
n = 1.5 in Fig. [3] shows the importance of the V term 
in the vertex Eq. ([4]). The non-interacting susceptibil- 
ity has a deep minimum on the Brillouin zone boundary. 
The charge response function on the other hand shows 
two different maxima at wave- vectors Qi and Q2. The 
CDW modulation related to these wave vectors are il- 
lustrated in Fig. [5] We will see that these instabilities 
occur over a wide area in the U — V plane. In fact apart 
from the phase separation instability (q = 0), which oc- 
curs for negative V , they are the only charge instabilities 
that can be found at this density. The situation changes 
as we change the density and one can expect to find in- 
commensurate CDW instabilities in another region of the 
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FIG. 6: (Color online) The charge response function at n = 
1.5 and T — 0.4 for different value of U and V. "When not 
specified, U takes the value U — 4. 
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FIG. 7: (Color online) The spin response function at fixed 
n = 1.5, f/ = 4, r = 0.4 and different values of F = 
0, 1, 1.5 and 1.8. 



U — V plane. 

Moving on to the spin susceptibility. Fig. [7] shows that 
the presence of the V term suppresses the spin response 
function more and more as V increases, concomitant with 
the increase in the charge response function. In princi- 
ple, one cannot find a strong maximum in both the spin 
and charge response functions. This is true in all one 
band homogeneous paramagnetic systems as dictated by 
Eqs. dll), dll) and g^^(O) = (Pauh sum-rule). Indeed, 
the Pauli sum-rule (obtained from the Pauli principle 
< >=< ria >) connects (7cc(0) to 5ss(0) in such a way 
that an increase in one forces a decrease in the other—. 



IV. CROSSOVER DIAGRAMS 

In mean-field theory, one normally finds finite tempera- 
ture phase transitions, in contradiction with the Mermin- 
Wagner theorem. In ETPSC, we obtain instead at a tem- 
perature Tx below which the correlation length begins to 
grow exponentially, diverging only at zero temperature. 
Tx is lower than the mean-field transition temperature 
because of the quantum Kanamori-Briickner renormal- 
ization of the vertices. In this low temperature regime, 
the characteristic frequency of the growing fiuctuations 
becomes less than temperature in dimensionless units. 
This is the renormalized-classical regime. Either the spin 
or the charge correlations grow exponentially at some 
characteristic wave vector that suggests which long-range 
order will likely be stabilized at zero temperature. Since 
our approach is not valid deep in the renormalized classi- 
cal regime, one cannot be sure that the zero-temperature 
phase will be precisely that suggested by the behavior at 
Tx- 

The value of Tx depends on density, U and V. We use 
x{<lx,<iv,0)/xo{cLx,q.y,0) = const to estimate Tx- For 
the sake of computational efficiency, we chose the con- 
stant to be 10 and checked that the general features do 
not change by choosing a larger value. This occurs be- 
cause the exponential growth of the correlations is rather 
sudden. A detailed discussion of this issue can be found 
in previous publications^^'^^. 

We present our results for Tx in Figs. [5] to [T3] as color 
(grey scale) plots in various planes of parameter space. 
There are four plots that present the U ~ V dependence 
of Tx at four densities, then two plots for the V — n 
dependence at fixed U. We indicate by lines of various 
colors and types the boundaries between regions where 
there is either a change in the wave vector of the growing 
correlations, or a change in the type of correlation, spin or 
charge. When we indicate a paramagnetic (Fermi liquid) 
region (PM), we mean that correlations did not grow, in 
either the spin or charge sectors, at temperatures as low 
as T = 0.01. 

Fig. [5] displays Tx at n = 1.5 as a function of U and V 
for both positive and negative values. At negative values 
of either V or U, superconducting correlations will be 
competing. Since superconductivity has not been taken 
into account here, the results in all quadrants, except 
the first one, should be taken as just indicative of what 
may happen in the spin or charge sectors. When V is 
negative, unless U is large, there is a strong tendency 
to phase separation (PS), i.e. the static charge response 
function starts, at Tx, to grow exponentially for wave 
vector q = 0. At positive U and V, incommensurate 
spin density waves (ISDW) are dominant, but U and V 
must be large enough, as expected from the absence of 
perfect nesting. At small U and V, the system remains 
paramagnetic (PM). Charge density waves appear at pos- 
itive U and V only if V is relatively large. Recall however 
that the effect of V is amplified by the presence of sev- 
eral neighbors. The charge instability in this parameter 
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FIG. 8: (Color online) Value of the crossover temperature Tx 
to the renormalized classical regime as a function of U and V 
at filling n = 1.5. The wave vector and spin or charge char- 
acter of the growing correlations is indicated by initials: In- 
commensurate spin density waves (ISDW), phase separation, 
or q = charge instability (PS), and incommensurate charge 
density waves (ISDW). The wave vectors of the two special 
charge density waves CDWl and CDW2 (both a/(3) x ^3)) 
are shown in Figs. [T] and O The color scale (grey scale) ap- 
pears on the right of the plot. Regions where either U or V 
are negative are shown for illustrative purposes only. 



range is of the CDWl type illustrated in Fig. [51 Charge 
instabilities are further amplified at positive V only if U 
is allowed to become negative. The y^S) x y^S) CDW2 
pattern is allowed only in extreme conditions of large 
positive V and large negative U . This is not suprising 
since one can see from the non-interacting Fermi surface 
in Fig. [1] that the corresponding wave vector Q2 is not 
particularly favored by nesting. The CDW2 phase is re- 
ally governed by properties of the vertex V , not so much 
by the non-interacting Fermi surface. 

When density is decreased to n = 4/3, the non- 
interacting Fermi surface becomes almost circular so at 
positive U and V the tendency to order is strongly sup- 
pressed, as can be seen from Fig. [31 Compared with 
the previous figure, the CDW2 vertex related instability 
is more robust while the CDWl and ICDW instabilities 
occur in smaller regions, the ICDW existing over a larger 
region this time than CDWl. At larger fillings, n = 5/3, 
where again the Fermi surface becomes almost circular, 
similar features are observed. 

As the filling decreases, the Fermi surface becomes 
more and more circular. Restricting ourselves to positive 
U and V , and staying at weak to intermediate coupling 
where our theory is valid, nothing interesting occurs. The 
system remains paramagnetic down to T = 0.01. We thus 
also present, in Fig.fTUland Fig. [Til results at large U and 
V where our theory is not strictly controlled. We feel 
these results are nevertheless interesting for two reasons. 




FIG. 9: (Color online) Crossover temperature Tx to the 
renormalized classical regime as a function of U and V at 
filling n = 4/3. Other symbols are defined in the caption of 
Fig. [5] Regions where either U or V are negative are shown 
for illustrative purposes only. 



First, some of our "phase boundaries" compare favorably 
with results obtained from other methods. Second, the 
renormalized classical regime occurs at such high tem- 
perature that U /T and V/T may begin to control the 
approach. 

Fig. [TUl and Fig. [Til thus show the crossover diagram 
for, respectively, n — 2/3 and 0.5 over a wide range of 
positive U and V. The CDW2 region now appears at 
positive U and V , contrary to the results in the previous 
figures, as long as the stabilizing interaction V is large 
enough. The boundary that separates CDW2 from PM 
in Fig. [TUl is very close to QMC results^^, which gives us 
confidence in the validity of the results. The ferromag- 
netic phase is dominant when both U and V are large. 
The ICDW phase does not appear at filling n = 2/3 
(Fig. [TTl) . The CDW2 phase is influenced to some extent 
not only by the vertex, but also by commcnsurability, as 
can be seen from the fact that it is more important at 
n — 2/3 than at n — 0.5. There are competing tenden- 
cies for the CDW2 phase: 1) The density n = 2/3 is more 
favorable for the commensurate CDW2 as reflected in the 
free response function, 2) The effect of off-site interaction 
is less important at lower value of the density. Based on 
these observations, one might surmise the presence of an 
optimal density where CDW2 appears over a larger area 
of U — V space. 

To explore in more details the density dependence, we 
present results as a function of n and V at fixed U. 
This was studied in particular by Motrunich and Lee^^i^i. 
They calculated the phase diagram with different meth- 
ods: a) renormalized mean field theory and variational 
quantum Monte Carlo with a trial wave function^^ and 
b) slave boson mean field theory^!. They suggest that 
the effect of the V term is taken into the account more 
accurately in the first method than in the second one. 
In the first paper, they found that the CDW2 phase can 
be reached at smaller V at the densities n = 1/3 and 
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FIG. 10: (Color online) Crossover temperature Tx to the 
renormalized classical regime as a function of U and V at 
filling n — 2/3. Other symbols are defined in the caption 
of Fig. |5] The only new acronym, FM stands for ferromag- 
netic. Results for the regions where U and AV are larger than 
about half the bandwidth are presented only for illustrative 
purposes. 



FIG. 12: (Color online) Crossover temperature Tx to the 
renormalized classical regime as a function of V and n at U = 
5. The yellow part of the figure represents a paramagnetic, 
or Fermi liquid, regime. Figs. [S] and [TO] for the meaning of 
symbols. 




FIG. 11: (Color online) Crossover temperature Tx to the 
renormalized classical regime as a function of U and V at 
filling n — 0.5. Other symbols are defined in the caption of 
Figs. |8]and[T0l Results for the regions where U and 41/ are 
larger than about half the bandwidth are presented only for 
illustrative purposes. 



2/3 than at other densities. This is not the case in the 
second paper where these densities play no special role. 
Since their calculations are at J7 ^ oo, this suggests that 
the general features of the phase diagram can be under- 
stood physically at low and high density: One needs a 
large V to stabilize the CDW at low value of the den- 
sity. Since at large value of the density, the hole plays 
the same role as the electron at low density, we expect 
at large V to find the CDW as well. This argument is 
correct if other phases do not suppress the CDW. 

We present the results of calculations at finite values 
of U in Figs. [12] and [HI In Fig. [HI the results in the 
V ~ n plane are for [/ = 5. It can easily be seen that 
CDW2 appears when both V and n are large. A ICDW 
region separates CDW2 from the paramagnetic, or Fermi 



liquid phase. There is a small area on the left of the 
plane where the ferromagnetic phase is stable. Fig. [T51 
shows the results for U = 10. CDW2 is still stable in the 
same region of the V — n plane but the small FM region 
of the previous figure has now grown and pushed away 
slightly the CDW2 phase. In other words, at larger U, 
spin fluctuations are playing a more important role, as 
expected. This is also shown by the appearance of an 
ISDW regime. It is obvious from Figs. [T^] and [T3] that 
the densities n = 1/3 and 2/3 do not play any special 
role, at least at these values of U. This is in agreement 
with the results of Ref. i24{. However, at large U, we 
find more ferromagnetic spin fiuctuations, a possibility 
that was not considered in Ref. [l^- We also calculated 
the same phase diagram at higher value of U where we 
find that the CDW2 region is completely swept away by 
ferromagnetism. 



V. FERMI SURFACE AND PSEUDOGAP 

It has already been shown that spin and superconduct- 
ing thermal fluctuations can2£i2ii2^iS, in two dimensions, 
open up a pseudogap on the Fermi surface that reflects 
the wave vector of the fluctuations. The same study can 
be performed here to check the effect of charge density 
wave fluctuations. We obtain the self-energy by substi- 
tuting our results for the susceptibilities in Eq. (fTT|) , from 
which one can compute the spectral function. 

We begin by the effect of spin fluctuations. In Fig. [Til 
we plot A{q,uj = 0) for n = 1.5, U ^ 4, V = and 
T = 0.4. In the jargon, this is known as a Momentum 
distribution curve (MDC). The dashed Green line is the 
Brillouin zone. Following the largest intensity regions, 
one can recognize the shape of the non-interacting Fermi 
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FIG. 13: (Color online) Same as Figs. [12] but for V = 10. 



FIG. 15: (Color online) The spectral function as a function of 
Lo at n = 1.5, (7 = 4, V = and T = 0.4 for different values 
of the wave-vector. Symmetry points are defined in Fig. \T\ 
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FIG. 14: (Color online) The spectral function in the Fermi 
liquid regime as a function of wave-vector at uj — o, n = 1.5, 
[7 = 4, y = and T = 0.4. The dashed Green line is the 
Brillouin zone. 



surface illustrated for n = 1.5 in Fig. [TJ It is clear from 
this figure that the effect of the on-site interaction at this 
temperature is just to introduce damping. Since this is a 
region where ISDW appear at low temperature, the Fermi 
surface can be destroyed by lowering the temperature or 
increasing U. 

For better understanding, we plot in Fig. [15] the spec- 
tral function as a function of frequency at different wave 
vectors for the same parameters as in the previous fig- 
ure. These are Energy dispersion curves (EDC). One 
can clearly observe the quasi-particle dispersion relation, 
the effect of U appearing as damping. The extra features 
at higher frequency are precursors of extra bands that 
would appear in the ordered state. Indeed, the EDC are 
plotted in a regime where the correlation length asso- 
ciated with incommensurate fluctuations is three to four 
lattice spacings {x{lx,q.v,0)/xo{q.x,q.y,0) = 10), enough 
to enter the renormalized classical regime. A simple gen- 



eralization of an argument presented earlier— shows that 
in that regime, the self-energy can be approximated by 



A2 



iT 



(13) 



where the sum runs over all the equivalent maxima of 
the susceptibility, six of them for the present case. Sub- 
stituting this expression into the general result for the 
spectral weight 



-2S" (k,w) 



(c^-ek-S' (k,c^))" + E" (k,c^) 



(14) 



one can reproduce qualitatively the behavior in Fig. 1151 
With larger correlation length (larger A), the spectrum 
clearly splits into more bands. 

The spin fluctuation induced precursor effects are also 
illustrated in Fig. [TC] that displays the density of states 
for various values of U and for densities near n = 1.5. 
The U = 2 curve is close to the non-interacting result. 
For stronger interaction, the extra band of states can be 
understood as a precursor effect by drawing the density 
of states for a static modulation. 

The fact that renormalized classical spin fluctuations 
can create a pseudogap on the Fermi surface in two di- 
mensions is well documented so we do not show a fig- 
ure corresponding to this case. We present the case of 
a Fermi surface pseudogap that originates from charge 
fluctuations. Increasing V enough that CDWl fluctua- 
tions become important, one can observe a pseudogap, 
as seen in A{q, uj — 0) for n = 1.5, U = 4, V = 1.95 and 
T — 0.4 displayed in Fig.[T7] The pseudogap opens up at 
the M point, in other words the spectral function at this 
point is suppressed. By analogy with the spin-fluctuation 
induced pseudogap, the symmetry equivalent M points 
are nearly connected by the symmetry equivalent wave 
vectors Ql of the CDWl charge instability. 
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FIG. 16: (Color online) The density of state for different val- 
ues of the densities and U for V = G, T = 0.4. 
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FIG. 17: (Color online) The spectral function as a function of 
wave-vector a.t uo = ^, n = 1.5, U — 4, V = 1.95 and T = 0.4 



The EDC corresponding to this situation are illus- 
trated in Fig. [181 It is clear again from this figure that 
the pseudogap opens up around the M point where the 
suppression occurs in the MDC of Fig. [iTl For other 
points the peak positions change only slightly. 



VI. DISCUSSION AND CONCLUSION 

We have used ETPSC to clarify the leading instabil- 
ities of the extended Hubbard model on the triangular 
lattice. The method is valid from weak to intermediate 
coupling. It satisfies the Mermin- Wagner theorem, in- 
cludes quantum (Kanamori-Briiuckner) renormalization 
of the vertices and does not assume a Migdal theorem for 
the self-energy. We interpret the entry into the renormal- 
ized classical regime (exponential growth of the correla- 



FIG. 18: (Color online) The spectral function as a function of 
Lo a.tn = 1.5, U = 4:,V = 1.95 and T = 0.4 for different value 
of the wave- vector. Symmetry points are defined in Fig. [1] 



tion length) as an indication of the phase that acquires 
long-range order at zero temperature. Since we scan all 
wave vectors for both spin and charge instabilities, our 
method is not biased towards a restricted set of instabili- 
ties as most other studies. Superconducting instabilities, 
however, have not been explored. 

The range of possible phases as a function of on-site 
interaction [/, nearest-neighbor interaction V and filling 
n is quite rich. Notwithstanding superconducting insta- 
bilities, negative values of U and V favor charge instabil- 
ities either at zero wave vector (phase separation) when 
U dominates, or at the CDWl and CDW2 wave vec- 
tors when V dominates. In the physically more relevant 
regime where U and V are both positive, spin instabilities 
are favored when U dominates and they occur at wave 
vectors that are essentially determined by the Fermi sur- 
face. That is particularly clear at large fillings where the 
shape of the Fermi surface is non-trivial. Larger V, on the 
other hand, favors charge instabilities that are generally 
determined by the vertex itself instead of by details of the 
Fermi surface. This is particularly clear at small filling 
where the Fermi surface is essentially circular. That pre- 
dominance of the vertex is why the CDW2 ( v^3) x -\/(3)) 
phase, for example, exists over a wide range of fillings 
and is not favored by commensurate fillings. We find 
that this CDW2 phase can compete with spin instabili- 
ties, especially ferromagnetism, when U is large as well. 
Ferromagnetic fiuctuations appear in a range of doping 
similar to that observed for the cobaltates. That com- 
petition with ferromagnetism has not been taken into 
account in earlier studies^^i^i. All the phases, except the 
Fermi liquid one, appear at large values of the interac- 
tions when filling is small, somewhat outside the regime 
of validity of our approach. Nevertheless, agreement with 
other approaches suggest that ETPSC extrapolates in a 
reasonable way towards strong coupling. The disappear- 
ance of all phases except the Fermi liquid one at small 
coupling and small fillings (n = 2/3) on the triangular 
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lattice is a clear manifestation of the effects of frustra- 
tion, as can be seen by contrasting with the square lattice 
case^^ {n = 0.75) where this does not occur. 

Finally, we also showed that charge-density wave ther- 
mal fluctuations can also induce a pseudogap. A pseudo- 
gap associated with CDW is observed2.S. experimentally. 
If its origin is the one discussed in the present paper, it 
should disappear as temperature rises above that where 
the charge correlation length becomes of the order of the 
thermal de Broglie wave length^, as observed in the spin 
fluctuation case in electron-doped cuprates^. 

Further studies should include the competition with 
superconducting fluctuations as well as more realistic 
modeling of the specific materials to which one wishes 
to apply our results. For example, the n — 1 case is 
relevant for the layered organics. 
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